NON-NEGATIVE INTEGER LINEAR CONGRUENCES 
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Abstract. We consider the problem of describing all non- negative integer solutions 
to a linear congruence in many variables. This question may be reduced to solving the 
congruence X\ + 2x 2 + 3x 3 + . . .+ (n — l)x n _x = (mod n) where eN = {0, 1,2,.. .}. 
We consider the monoid of solutions of this equation and prove a conjecture of Elashvili 
concerning the structure of these solutions. This yields a simple algorithm for generating 
most (conjecturally all) of the high degree indecomposable solutions of the equation. 



1. Introduction 

We consider the problem of finding all non-negative integer solutions to a linear con- 
gruence 

W\Xi + w 2 x 2 + . . . + w r x r = (mod n) 

By a non-negative integer solution, we mean a solution A = (ai,a 2 , . . . , a r ) with G 
N:= {0,1,2,...} for all i = 1, 2, . . . , r. 

As one would expect from such a basic question, this problem has a rich history. The 
earliest published discussion of this problem known to the authors was by Carl W. 
Strom in 1931 ( jStlj ). A number of mathematicians have considered this problem. 
Notably Paul Erdos, Jacques Dixmier, Jean-Paul Nicolas ( jDENj ). Victor 
Kac, Richard Stanley ( K ) and Alexander Elashvili ([Ej). 

In particular, Elashvili performed a number of computer experiments and made a 
number of conjectures concerning the structure of the monoid of solutions. Here we prove 
correct one of Elashvili's conjectures. This allows us to construct most (conjecturally 
all) of the "large" indecomposable solutions by a very simple algorithm. 

Also of interest are the papers jEJlj . |EJ2j by Elashvili and Jibladze and [EJP 
by Elashvili, Jibladze and Pataraia where the "Hermite reciprocity" exhibited by 
the monoid of solutions is examined. 

2. Preliminaries 

We take N = {0, 1,2,.. .} and let n be a positive integer. Consider the linear congru- 
ence 

(2.0.1) WiXi + w 2 x 2 + . . . + w r x r = (mod n) 

where Wi, w 2} ■ ■ ■ , w r G Z and unknowns. We want to describe all 

solutions A = (ai, a 2 , ... , a n -i) G N r to this congruence. 
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Clearly all that matters here is the residue class of the Wi modulo n and thus we may 
assume that < Wi < n for all i. Also if one of the W{ is divisible by n then the equation 
imposes no restriction whatsoever on Xj and thus we will assume that 1 < < n for all 
i. 

If w\ = w 2 then we may replace the single equation (j2.0.1J) by the pair of equations 

Wiyi + W3X3 + . . . + w r x r = (mod n) and X\ + x 2 = y\. 

Thus we may assume that the Wi are distinct and so we have reduced to the case where 
{w\, . . . , w r } is a subset of {1, 2, n — 1}. Now we consider 

(2.0.2) xi + 2x 2 + 3x 3 + . . . + (n - l)x n _i = (mod n) 

The solutions to (|2.(J.1|) are the solutions to (|2.(J.2|) with X{ = for all i ^ {wi, . . . , w r }. 
Hence to solve our original problem it suffices to find all solutions to Equation ()2.().2|) . 

3. Monoid of Solutions 

We let M denote the set of all solutions to Equation ()2.0.2|) . 

M := {x G N"^ 1 I x x + 2x 2 + . . . + (n - l)z n -i = (mod n)} . 

Clearly M forms a monoid under componentwise addition, i.e., M is closed under this 
addition and contains an additive identity, the trivial solution = (0, 0, . . . , 0). 

In order to describe all solutions of ()2.0.2|) explicitly we want to find the set of minimal 
generators of the monoid M. We denote this set of generators by IM. We say that a 
non-trivial solution A G M is decomposable if A can be written as non-trivial sum of two 
other solutions: A = B + C where B, C 7^ 0. Otherwise we say that A is indecomposable 
(also called non-shortenable in the literature). Thus IM is the set of indecomposable 
solutions. 

We define the degree (also called the height in the literature) of a solution A = 
(ai,a 2 , . . . ,a n _i) G M by deg(A) = a\ + a 2 + . . . + a n ~\ and we denote the set of 
solutions of degree k by M(k) := {A G M \ deg(A) = k}. Similarly, we let IM{k) 
denote the set of indecomposable solutions of degree k: IM{k) = IM R M{k). 

Gordan's Lemma [Gj] states that there are only finitely many indecomposable solutions, 
i.e., that IM is finite. This is also easy to see directly as follows. The extremal solutions 
Ei := (n, 0, ...,0), E 2 := (0, n, 0, . . . , 0), . . ., E n _i := (0, 0, . . . , 0, n) show that any 
indecomposable solution, (ai, a 2 , ... , a n ) must satisfy <n for all i. 

In fact, Emmy Noether [N] showed that if A is indecomposable then deg(v4) < n. 
Furthermore A is indecomposable with deg(A) = n if and only if A is an extremal 
solution Ei with gcd(i,n) = 1. For a simple proof of these results see jS]- 

We define the multiplicity of a solution A, denoted m(A) by 

. ai + 2a 2 + . . . + (n - l)a n -i 

m(A) := . 

n 

Example 3.1. Consider n = 4. Here /M = = (4,0,0),A 2 = (0,2,0),A 3 = 

(0,0, 4), At = (1,0, 1), A = (2, 1,0), A 6 = (0,1,2)}. The degrees of these solutions 
are 4, 2, 4, 2, 3, 3 respectively and the multiplicities are 1, 1, 3, 1, 1, 2 respectively. 
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4. The Automorphism Group 

Let G := AutiTLjriE). The order of G is given by <p(n) where <ft is the Euler phi 
function, also called the totient function. The elements of G may be represented by 
the 4>(n) positive integers less than n and relatively prime to n. Each such integer g 
induces a permutation, o = cr g , of {1, 2, . . . , n — 1} given by a{i) = gi (mod n). Let 
A = (oi, a 2 , . . . , a„_i) G M, i.e., a x + 2a 2 + . . . + (n — l)a n -i = (mod n). Multiplying 
this equation by g gives (g)ai + (2g)a 2 + (3g)a 3 + . . . + (gn — g)a n ^i = (mod n). 
Reducing these new coefficients modulo n and reordering this becomes a a -i^ + 2a a -i{2) + 
. . . + (n — l)ao-i( n _i) = (mod n). Thus if A = (ai, a 2 , . . . , a n _i) G M then g • A :— 
(a CT -i(i), 0,7-1(2), • • • , a CT -i( n _i)) G M 

Since g • A is a permutation of A, the action of G on M preserves degree, and thus 
G also acts on each M{k) for 6 N. Note however that the action does not preserve 
multiplicities in general. Furthermore if g G G and A = B + C is a decomposable 
solution, then g-A = g- B + g- C and therefore G preserves JM and each IM(k). 

Example 4.1. Consider n = 9. Here G is represented {1, 2, 4, 5, 7, 8} and the correspond- 
ing six permutations of Z/9Z are given by 0\ = e, a 2 = (1, 2, 4, 8, 7, 5)(3, 6), 04 = a\ = 
(1,4,7)(2,8,5)(3)(6), a 5 = a\ = (1, 5, 7, 8, 4, 2)(3, 6), a 7 = a\ = (1, 7, 4)(2, 5, 8)(3)(6) 
and cr 8 = cr| = (1, 8)(2, 7)(3, 6), (4, 5). Thus, for example, 2- (ai, a 2 , 03, a 4 , a 5 , a 6 , a 7 , a 8 ) = 
( a 5; a i; a 6! a 25 a 7) a 3> a 8> °-4) and 4-(ai, a 2 , 03, a.4, 05, a6, 07, as) = (07, 05, (33, a%, cig, 0.6, 04, a 2 ). 

Note that G always contains the element n — 1 which is of order 2 and which we 
also denote by —1. This element induces the permutation er_i which acts via —1 • 
(ai, a 2 , . . . , a n _i) = (a n _i, a n _ 2 , . . . , a 3 , a 2 , ai). 

Let F(n) denote the number of indecomposable solutions to Equation (|2.(J.2|) . F(n) : = 
#/M. Victor Kac jK] showed that the number of minimal generators for the ring of 
invariants of SL(2, C) acting on the space of binary forms of degree d exceeds F(d — 2) 
if d is odd. Kac credits Richard Stanley for observing that if A is a solution of 
multiplicity 1 then A is indecomposable. This follows from the fact that the multiplicity 
function m is a homomorphism of monoids from M to N and 1 is indecomposable in 
N. Kac also observed that the extremal solutions E$ (defined in Section above) with 
gcd(z,n) = 1 are also indecomposable. This gave Kac the lower bound F(n) > p(n) + 
4>(n) — 1 where p(n) denotes the number of partitions of n. 

Jacques Dixmier, Paul Erdos and Jean-Louis Nicholas studied the function 
F(n) and significantly improved Kac's lower bound QDENJ). They were able to prove 
that 

n l/2 

-p(n) 



lim inf F(n) 



> 



log n ■ log log n 

It is tempting to think that the 67-orbits of the multiplicity 1 solutions would comprise 
all elements of IM. This is not true however. Consider n = 6. Then G is a group of 
order 2, G = {1, -1}. The solutions A x = (1, 0, 1, 2, 0) and A 2 = -1 • A x = (0, 2, 1, 0, 1) 
are both indecomposable and both have multiplicity 2. 

We define the level of a solution A, denoted £(A), by £(A) = mm{m(g(A)) \ g G G}. 

Notethatm(A)+m(-l-A) = deg(A). This implies 2 YIbgg-a m {B) = deg(A)#(G ■ A) , 
i.e., that the average multiplicity of the elements in the G-orbit of A is half the degree 
of A. 
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5. Elashvili's conjectures 



In jEj, Elashvili made a number of interesting and deep conjectures concerning the 
structure of the solutions to Equation (|2.0.2j) . In order to state some of these conjectures 
we will denote by p(t) the number of partitions of the integer t. We also use \n/2\ to 
denote the greatest integer less than or equal to n/2 and [n/2] := n — [n/2\. 

Conjecture 1: If A G IM{k) where k > [n/2\ + 2 then £(A) = 1. 

Conjecture 2: If k > \n/2\ + 2 then IM(k) consists of p(n — k) orbits under G. 

Conjecture 3: If k > \n/2\ + 2 then IM(k) contains exactly p(n — k) orbits of level 1. 

Here we prove Conjecture 3. Furthermore we will show that if k > [n/2] + 1 then 
every orbit of level 1 contains exactly one multiplicity 1 element and has size <p(n). Thus 
if k > [n/2] +1 then IM(k) contains exactly (p(n)p(n — k) level 1 solutions. 

This gives a very simple and fast algorithm to generate all the level 1 solutions whose 
degree, k, is at least [n/2] + 1 as follows. For each partition, n — k = b\ + b 2 + • • • + b s , 
of n — k put b s+ i = ■ ■ ■ = bk = and define q := b,i + 1 for 1 < % < k. Then define A via 
a i '■= #{j : c j — 0- This constructs all multiplicity 1 solutions if k > [n/2] + 1. Now 
use the action of G to generate the </>(n) solutions in the orbit of each such multiplicity 
1 solution. 

If Conjecture 2 is true then this algorithm rapidly produces all elements of IM(k) for 
k > [n/2\ +2. This is surprising, since without Conjecture 2, the computations required 
to generate the elements of IM(k) become increasingly hard as k increases. 



Before proceeding further we want to make a change of variables. Suppose then that 
A G M(k). We interpret the solution A as a partition of the integer m(A)n into k parts. 
This partition consists of a\ l's, a 2 2's,. . . , and a n _i (n-l)'s. We write this partition as 
an unordered sequence (or multi-set) of k numbers: 



The integers yi,y2,---,Vk with 1 < y^ < n — 1 for 1 < i < k are our new variables for 
describing A. Given [yi, y 2 , . . . , yk] we may easily recover A since := #{j | yj = i}. 
We have yi + y 2 + ■ ■ ■ + yk = m(A)n. 

Notice that the sequence y\ — 1, y 2 — 1, . . . , yk — 1 is a partition of m(A)n — k. Fur- 
thermore, every partition of m(A)n — k arises from a partition of m(A)n into k parts in 
this manner. 

The principal advantage of this new description for elements of M is that it makes the 
action of G on M more tractable. To see this let g G G be a positive integer less than n 
and relatively prime to n. Then g ■ [yx,y 2 , . . . ,y k ] = [gy x (mod n), gy 2 (mod n),...,gy k 
(mod n)]. 

Now we proceed to give our proof of Elashvili's Conjecture 3. 

Proposition 6.1. Let A G M(k) and let 1 < g < n — 1 where g is relatively prime to n 
represent an element ofG. Write B = g-A, andu = m(A) andv = m(B). If k > gu — v 
then ug 2 — (k + u + v)g + v(n + 1) > 0. 



6. Proof of Conjecture 3 
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Proof. Write A = [yi, y 2 , . . . , y k ] where y\ > y 2 > . . . > y k . For each i with 1 < i < k we 
use the division algorithm to write gyi = q^n + r« where qi G N and < r« < n. Then 
B = [r±, r 2 , ■ ■ ■ , r k ]. Note that the r\ may fail to be in decreasing order and also that no 
r\ can equal 0. 

Now gun = gfa + y 2 + . . . + y k ) = (qin + n) + (q 2 n + r 2 ) + . . . + (q k n + r k ) = 
(?i + ?2 + • • • + Qk)n + (ri + r 2 + . . . + r k ) where r 1 + r 2 + . . . + r k = vn. 
Therefore, gu = (qi + q 2 + . . . + q k ) + v. 

Since yi > y 2 > ■ ■ ■ > y k , we have q± > q 2 > ■ ■ ■ > q k . Therefore from gu — v = Yli=i Qi 
we conclude that ^ = for all % > gu — v. Therefore 



gu—v gu—v 



i=i 

ug—v 

9^2(Vi- l)+g(gu-v) 



i=i i=i 



i=l 
fc 



i=l 



= g(un — k) + g 2 w — gt> 



Also 



gu—v gu—v 



i=l 



1 i=l 

i=l 

> gwn — vn + gu — v 



Combining these formulae we obtain the desired quadratic condition ug 2 — (k + u + 
w)^ + w(n+ 1) > 0. □ 

Now we specialize to the case u — v — 1. Thus we are considering a pair of solutions 
A and B = g ■ A both of degree A; and both of multiplicity I. 

Lemma 6.2. Let A e M(/c) 6e a solution of multiplicity 1. Write A = [yi,y 2 , . . . ,y k ] 
where y 1 > y 2 > ■ ■ ■ > y k . If k > [n/2\ + 2 then y k _ 2 = y k ^ = y k = l. If k > \n/2] + 1 
then y k -i = y k = l. 

Proof. First suppose that k > [n/2\ + 2 and assume, by way of contradiction, that 
y k _ 2 > 2. Then n = (y x +y 2 + . . .+y k - 2 )+y k -i+y k > 2(fc-2) + l + l > 2[n/2j +2 > n + 1. 

Similarly if k > [n/2] + 1 we assume, by way of contradiction, that y k -\ > 2. Then 
n = (y 1 + y 2 + ... + y fe _ 1 )+y fc >2(A;-l) + l>2([n/2l) + l>n + l. ' □ 

Proposition 6.3. Let A e M{k) be a solution of multiplicity 1 where k > [n/2] + 1. 
TTien t/ie G-orbit of A contains no other element of multiplicity 1. Furthermore, G acts 
faithfully on the orbit of A and thus this orbit contains exactly <f>(n) elements. 
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Proof. Let B = g ■ A for some g G G and suppose B also has multiplicity 1. Lemma 16. 21 
implies that B = g ■ A = [r±, T2, ■ ■ ■ , r/-_2, <?, <?]■ Since I? has multiplicity 1, we have 
n = ri + r 2 + . . . + r^_2 + g + g>2g + k — 2 and thus g < (n — k + 2)/2 < k/2. From 
this we see that the hypothesis k > gu — v is satisfied. Therefore by Proposition 16. 1\ g 
and k must satisfy the quadratic condition 

g 2 -(k + 2)g + (n + l)>0 . 

Let / denote the real valued function f(g) = g 2 — {k + 2)g + (n + 1). Then /(l) = 
n — k > and /(2) = n+ l — 2/c<0 and thus / has a root in the interval [1,2). Since 
the sum of the two roots of / is k + 2 we see that the other root of / lies in the interval 
(k, k + 1]. Thus our quadratic condition implies that either g < 1 or else g > k + 1. But 
we have already seen that g < k/2 and thus we must have g = 1 and so A = -B. 

This shows that the G-orbit of A contains no other element of multiplicity 1. Further- 
more, G acts faithfully on this orbit and thus it contains exactly 0(n) elements. □ 

Remark 6.4. Of course the quadratic condition ug 2 — (k + u + v)g + v(n + 1) > can 
be applied to cases other than it = f = 1. For example, taking u = v = 2 one can show 
that a solution of degree k (and level 2) with > (2n + 8)/3 must have an orbit of size 
4>(n) or 0(n)/2. 
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